% Main_Framework of searching multiple solution points
% A Piecewise Linear dynamical model is investigated herein to detect multiple solution points;
% The Piecewise Linear model is symmetrical about the origin and is determined by K1_Array, K2_Array, B1_Array and B2_Array.
    tic;
	K1_Array = 5;
	K2_Array = 30; 
	B1_Array = 0.3;
	B2_Array = 0.5;
	KL_Array = 0;
	CL_Array = 0.1;
	gama_Array = 0.4;
	ws=1.5;
	NH=1;
	Piecewise_Data = [K1_Array,K2_Array,KL_Array,CL_Array,B1_Array,B2_Array,gama_Array];
	load 'Initial_Solution_Set.mat';
    Relative_Amplitude=zeros(size(Initial_Solution_Set,1),1);
    R_Amplitude=zeros(size(Initial_Solution_Set,1),length(ws));
	figure(1);hold on;
    for jndex = 1 : length(ws)
        w = ws(jndex);
        parfor index = 1 : size(Initial_Solution_Set,1)
            [Relative_Amplitude(index),~]=HB_Verification(w,NH,Piecewise_Data,Initial_Solution_Set',index);
        end
        R_Amplitude(:,jndex) = Relative_Amplitude;     
    end
    plot(ws,R_Amplitude','r*');         
    toc;